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The coefficients of the volume, surface, coulomb, asymmetry and pair- 
ing energy terms of the semiempirical liquid drop model mass formula 
have been determined by furnishing best fit to the observed mass excesses. 
Slightly different sets of the weighting parameters for liquid drop model 
mass formula have been obtained from minimizations of chisquare and mean 
square deviation. The most recent experimental and estimated mass ex- 
cesses from Audi-Wapstra-Thibault atomic mass table have been used for 
the least square fitting procedure. Equation of state, nuclear incompress- 
ibility, nuclear mean free path and the most stable nuclei for corresponding 
atomic numbers, all are in good agreement with the experimental results. 
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1. Introduction 

A mass formula is not only useful to the evaluation of atomic masses and 
nuclear binding energies but provides theoretical predictions concerning a 
number of features of nuclei and their behaviour. Mass models started out 
as empirical formulas and have now evolved into more or less full scale nu- 
clear models, many capable of predicting a whole host of nuclear properties 
besides the mass. The atomic or nuclear mass is one of the most decisive 
factors governing nuclear stability. Enjoying reasonable success and relative 
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ease of use are the class of mesoscopic models, where a gross macroscopic 
component is simply dependent on the number of nucleons and is finely 
sculpted using microscopic corrections for shell behaviour, deformation, nu- 
clear incompressibility and nucleon pairing. This approach finds its orgin 
in the famous Bethe [Q-Weizsacker [2] mass formula which is an empirically 
refined form of the liquid drop model for the binding energy of nuclei. Since 
Bethe- Weizsacker mass formula is based on a liquid drop description of the 
nucleus, it is expected to reproduce the gross features of nuclear binding 
energies better for medium and heavy nuclei than for light nuclei. 

Different approaches to the evaluation of the coefficients of the volume, 
surface, coulomb, asymmetry and pairing energy terms of the semiempiri- 
cal Bethe- Weizsacker mass formula furnishing the best fit to the observed 
masses, have yielded different sets of results. Some of which are compared 
in Refs. |H1[3J- The mass of an atom and especially of a nucleus is not merely 
the sum of the masses of its constituents but is accompanied by a binding 
energy. In nuclear physics accurate mass measurements have a great influ- 
ence as the binding energy reveals a host of nuclear structure effects, magic 
configuration number effects known as shell and subshell closures, shape 
and deformation effects and also nucleon pairing. Finally, as the binding 
energy determines how much energy is available for a given nuclear reaction, 
the impact of masses in nuclear astrophysics is also far reaching. A host of 
quantities derived from the masses are in play, notably, neutron and pro- 
ton separation energies. This atomic mass evaluation currently containing 
almost harmoniously adjusted values is published every few years and the 
recent one is available in reference [3] . In the present work, a five parameter 
least square fit to the Bethe- Weizsacker mass formula has been achieved for 
new evaluations of energy coefficients by minimizing the chisquare and the 
mean square deviation for atomic mass excesses using the latest compilation 
[3] and its consequences on nuclear properties such as incompressibility of 
nuclear matter, Coulomb radius constant and effect of asymmetry energy 
coefficient have been investigated. 

2. The liquid drop model nuclear binding energy and the atomic 

mass excess 

The Bethe- Weizsacker formula H] is an empirically refined form of 
the liquid drop model for the binding energy of nuclei. Expressed in terms 
of the mass number A and the atomic number Z for a nucleus, the binding 
energy B(A,Z) is given by, 



B(A, Z) = a v A - a s A 2/3 - a c Z(Z - l)/A 1/3 - a sym (A - 2Z) 2 /A + 5, (1) 
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where 5 = a p A ' for even N — even Z, 
= -a p A~ 1/2 for odd N - odd Z, 
= for odd A, (2) 

and N is the neutron number of the nucleus. 

The above formula for nuclear binding energy was prescribed by Bethe 
and Weizsacker [2j where the nucleus was considered as a droplet of 
incompressible matter and all nuclei have the same density. Scattering ex- 
periments suggest that nuclei have approximately constant density, so that 
the nuclear radius can be calculated by using that density as if the nucleus 
were a drop of a uniform liquid. Thus the first approximation to the bind- 
ing energy was identified as due to the saturated exchange force and the 
term a v A was called the volume energy. The second term — a s A 2 / 3 , called 
the surface energy, was a correction to it since a deficit of binding energy 
for surface nucleons is expected due to those nucleons which are visualized 
as being at the nuclear surface have fewer near neighbours than which are 
deep within the nuclear matter. Only well known infinite range force in a 
nucleus is the Coulomb repulsion among protons that gives rise to the re- 
pulsive Coulomb energy term — a c Z {Z-\)IA X I Z . This term was derived by 
calculating the Coulomb energy assuming that the total nuclear charge Ze 
to be spread uniformly throughout the spherical nuclear volume of radius 
equal to r^A 1 ^ , ro being the Coulomb radius constant. The term Z(Z — 1) 
instead of Z 2 appears due to the fact that the nuclear charges are integral 
multiples of electronic charge e and the single proton should not have any 
contribution to coulomb energy and hence Coulomb self energy contribution 
by each of Z protons should be removed. The Coulomb energy coefficient 
a c is given by 



a c = (3e 2 /5r ). (3) 

Another deficit of binding energy depends on the neutron excess N — Z = 
A — 2Z since these excess neutrons have to go into previously unoccupied 
quantum states having larger kinetic energy and smaller potential energy 
than those already occupied. This asymmetry energy term — a sym ^ A ~^ , 
thus, arises due to purely quantum mechanical effect. The facts that there 
is a finite pairing energy between odd- A and even- A nuclei and the anoma- 
lously large binding energy for nuclei containing a magic number of neutrons 
and protons, fail to appear in liquid drop model since intrinsic spins of the 
nucleons have been omitted. To correct for this omission the pairing energy 
term 5 was empirically incorporated. 
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Using the definition of nuclear binding energy B(A,Z) which is defined 
as the energy required to separate all the nucleons consituting a nucleus, 
the mass equation can be expressed as 



where m p and m n are the masses of proton and neutron respectively and 
M nuc i eus (A, Z) is the actual mass of the nucleus. 

The reason that the atomic mass is considered rather than the nuclear 
mass is that historically, the former has been the actual experimentally 
measured quantity, whereas the latter is less accurate bacause its extraction 
requires a knowledge of binding energy of the Z atomic electrons. However, 
recent developments now allow the nuclear masses to be measured directly 
[Hj. For those applications where it is necessary to know the actual mass 
of the nucleus M nuc [ eus (A, Z) itself, its value (in MeV) can be found from 
Eq.(4) using the nuclear binding energy or from the atomic mass excess by 
use of the relationship 



M nudeus (A, Z)=Au + AMa,z ~ Zm e + a el Z 2 - 39 + b el Z™ 5 , (5) 



where m e is the mass of an electron, the atomic mass unit u is 1/12 the 
mass of 12 C atom, AMa,z is the atomic mass excess of an atom of mass 
number A and atomic number Z and the electronic binding energy constants 
[7j a e i = 1.44381 x 1(T 5 MeV and b d = 1.55468 x 1(T 12 MeV. Hence from 
Eq.(4) and Eq.(5) the atomic mass excess is given by 

AM A ,z = ZAm H + (A — Z)Am n - a d Z 2m - b el Z 5 - 35 - B(A, Z) (6) 

where Amu = m p +m e —u = 7.28897050 + a e i + b e i MeV and Am n = m n —u 
= 8.07131710 MeV. 

3. The mass fitting and the extraction of energy coefficients 

Different approaches to the evaluation of the weighting parameters a v , 
a s , a c , a sym and a p furnishing the best fit to the observed masses when 
inserted in the mass formula have yielded different sets of results. Some of 
which are compared in Refs. [31 H]. The two quantities of importance that 
can be used for least square fitting are the root mean square deviation a 
and x 2 /N which are defined as 



M nuc i eus (A, Z) = Zm p + (A - Z)m n - B(A, Z) 



(4) 




(7) 
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and X 2 /N= (l/JV)^[(AM n -AM Ex )/AM E J 2 (8) 

respectively where the summations extend to iV data points for which mea- 
sured atomic mass excesses and corresponding binding energies are known. 

Table 1. Energy coefficients for the Bethe-Weizsacker mass formula. 



a v 


a s 


a c 


O'sym 


a p 


MeV 


MeV 


MeV 


MeV 


MeV 


15.260 


16.267 


0.689 


22.209 


10.076 


±0.020 


±0.062 


±0.001 


±0.048 


±0.854 


[a] 


Minimized 


quantity 


a 2 = 


15.15 


15.777 


18.340 


0.710 


23.210 


11.998 


±0.053 


±0.174 


±0.003 


±0.103 


±1.912 


[a] 


Minimized 


quantity 


X 2 /N= 


3.336 


15.409 


16.873 


0.695 


22.435 


11.155 


±0.026 


±0.080 


±0.002 


±0.065 


±0.864 


[b] 


Minimized 


quantity 


a 2 = 


11.74 


15.777 


18.341 


0.710 


23.211 


11.996 


±0.037 


±0.133 


±0.002 


±0.060 


±1.536 


[b] 


Minimized 


quantity 


X 2 /N= 


2.484 



[a] . Using both measured 2228±extrapolated 951 atomic mass excesses. 

[b] . Using only the experimentally measured 2228 atomic mass excesses. 

The quantity AM^ is theoretically calculated atomic mass excess obtained 
using Eq.(6) while AMe x is the corresponding experimental atomic mass 
excess obtained from the Audi-Wapstra-Thibault atomic mass table [Sj. For 
the calculations of x 2 /N and the mean square deviation a 2 , masses p] 
of 3179 nuclei including the 951 extrapolated values which are predicted 
according to systematics, have been used for least square fitting. These 
calculations have been repeated using only the 2228 experimental data. The 
X 2 /iV and a 2 minimizations have been performed separately since both can 
not be minimized with identical sets of values for energy coefficients. In case 
of the x 2 /N minimization, the denominator has to be kept equal to unity 
only for 12 C in order to accommodate it in the calculation as the atomic 
mass excess of 12 C is, by definition, equal to zero. In table-1, values of the 
energy coefficients along with the minimized quantities with their values 
have been listed. Since Bethe-Weizsacker mass formula is more appropriate 
for nuclear binding energies of medium and heavy nuclei than for light nuclei, 



G 
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the values obtained for a 2 and x 2 /N are not very small. However, as can 
be seen from table- 1, the errors associated with the energy coefficients a v , 
a s , a c and a sym are very small while that with a p is reasonably small. 

As can be seen from the table- 1, both a 2 , which is a measure of ab- 
solute error, and the x 2 /N , which is a measure of relative error, can not 
be minimized simultaneously with the same set of parameters. The set of 
parameters obtained by minimizing x 2 /N are almost same as those listed in 
Refs. [HUH]. However, in many cases, such as reaction Q value calculations, 
absolute errors involved in atomic mass excesses are the quantities of concern 
and a 2 minimization plays more important role than x 2 /N minimization. 
On the contrary, x 2 /N minimization causes minimization of relative, that 
is, uniform percentage error involved in the mass predictions. 



4. Energy coefficients and the nuclear matter properties 

A density dependent M3Y effective nucleon-nucleon (NN) interaction 
based on the G-matrix elements of the Reid-Elliott NN potential has been 
used to determine the nuclear matter equation of state. The equilibrium 
density of the nuclear matter has been determined by minimizing the en- 
ergy per nucleon. The density dependence parameters have been chosen to 
reproduce the saturation energy per nucleon and the saturation density of 
spin and isospin symmetric cold infinite nuclear matter. The general ex- 
pression for the density dependent effective NN interaction potential v(s) is 
written as 

v(s,p,e) =t M3Y (s,e)g{p,e) MeV (9) 

where M3Y effective interaction potential supplemented by a zero range 
pseudopotential t M3Y is given by 

—4s —2.5s 

t M2Y (s, e) = 7999 (MeV). - A 2134 (MeV). - + Joo(e)S(s) (10) 

4s 2.5s 

where the costants 4 and 2.5 have dimensions of fm , 5(s) has the dimen- 
sion of fm~ 3 and the zero-range pseudo-potential Joo(e) representing the 
single-nucleon exchange term is given by 

J 00 ( e ) = -276(1 - ae)(MeV.fm 3 ) (11) 
and the dimensionless density dependent part is given by 



5 (p,e)=C(l-/3(e)p 2 / 3 ) 



(12) 



Nuclear printed on February 4, 2008 



7 



The energy per nucleoli e obtained using the effective nucleon-nucleon 
interaction v(s) for the spin and isospin symmetric cold infinite nuclear 
matter, which will henceforth be called the standard nuclear matter, is 
given by 

e = [3h 2 k 2 F /10m]+g(p,e)pJ v /2 = [3h 2 k 2 F /10m] + [pJ v C(l- f3p 2/3 )/2] (13) 

where m is the nucleonic mass, k F = (1.57T 2 p) 1 ^ 3 is the Fermi momentum, 
p is the nucleonic density while po being the saturation density for the 
standard nuclear matter and J v represents the volume integral of t M3Y , 
the M3Y interaction supplemented by the zero-range pseudopotential. The 
Eq.(13) can be differentiated with respect to p to yield equation 

Oe/dp = [h 2 k 2 F /5mp] + J V C[1 - (5/3)(3p 2/3 }/2 (14) 

The equilibrium density of the nuclear matter is determined from the satu- 
ration condition de/dp = 0. Then Eq.(13) and Eq.(14) with the saturation 
condition can be solved simultaneously for fixed values of the saturation 
energy per nucleon eo and the saturation density po of the standard nuclear 
matter, to obtain the values of the density dependence parameters @ and C 
given by 

/3=[(l-p) Po 2/3 ]/[3-(5/3)p], (15) 
p= [10me ]/[h 2 k Fo ], (16) 
k Fo = [1.5^ 2 po] 1/3 , (17) 

C = -[2h 2 k 2 Fo ]/[5mJ vPo (l - (5/3)pp 2 /3 )}, (18) 

respectively. It is quite obvious that the density dependence parameter (5 
obtained by this method depends only on the saturation energy per nu- 
cleon eo, the saturation density po but not on the parameters of the M3Y 
interaction while the other density dependence parameter C depends on the 
parameters of the M3Y interaction also through the volume integral J v . The 
energy per nucleon can be rewritten as 

6 = [3h 2 k 2 F /Wm] - {p/p Q )[h 2 k 2 Fo {\ - (3p 2 / 3 )]/[5m(l - (5/3)/3p^ /3 )] (19) 

and the incompressibility Kq of the spin and isospin symmetric cold infinite 
nuclear matter is given by 
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K = k 2 F d 2 e/dk 2 F = 9p 2 d 2 e/dp 2 \ p=po = [-(3h 2 k Fo /5m) - 5J v C/3pl /3 } (20) 
The pressure P and the energy density e of nuclear matter can be given by 
P = p 2 de/8p = [ph 2 k 2 F /bm] + p 2 J v C[l - (5/3)/3p 2 / 3 ]/2, (21) 

e = p{e + mc 2 ) = p[(3h 2 k 2 F /10m) + pJ v C(l - (3p 2/3 )/2 + mc 2 ], (22) 

respectively, and thus the velocity of sound v s in standard nuclear matter 
is given by 

= M1I>«^| < 23 > 

One of the most directly derivable information provided by the mass 
formula is the magnitude of the Coulomb radius constant contained within 
the expression for the Coulomb energy constant a c = (3e 2 /5ro) assuming 
a uniform volume distribution of charge. Using the values of Coulomb en- 
ergy coefficients a c , the values obtained for the Coulomb radius constant 
ro have been listed in table-2. The present value for the Coulomb radius 
constant is ro = 1.22 — 1.25/m which like previous such mass fittings still 
overpredicts. The crucial property of nuclear matter is the saturation den- 
sity. It has been recognised right from the beginning that Bethe-Weizsacker 
mass formula or its improved versions can not give this property through 
its Coulomb radius constant using the simple relation po = 3/47rrg. This 
density is measured through electron scattering experiments on heavy nu- 
clei. This gives a value corresponding to ro = 1.12 — 1.13 fm. Identifying a v 
as the saturation energy per nucleon for the spin and the isospin symmet- 
ric cold infinite nuclear matter and hence using the saturation energy per 
nucleon equal to — 15.260MeV along with the density dependent M3Y effec- 
tive interaction and the commonly used value for saturation density equal 
to 0.1533/m -3 |1U| . as one can see in table-2, the nuclear incompressibility 
is found to be 293AMeV which is in close agreement to experimental data 
[111 112] . In Fig. ^the energy per nucleon E/A=e of standard nuclear matter 
is plotted as a function of nucleonic density p. The continuous line repre- 
sents the curve using saturation energy per nucleon of —15.26MeV and the 
dots represent the same using saturation energy per nucleon of — 15.78MeV 
whereas the dash-dotted line represents the same for the A18 model using 
variational chain summation (VCS) ^5] f° r the spi n an d isospin symmetric 
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infinite nuclear matter. The minima of the energies per nucleon equaling the 
saturation energies of -15.260 MeV and -15.777 MeV for the present calcula- 
tions occur correctly at the saturation density po = 0.1533/m~ 3 while that 
for the A18(VCS) model occurs around p = 0.28/m~ 3 with a saturation 
energy of about -17.3 MeV. 

Table 2. Coulomb radius constant and the compression modulus using the energy 
dependence parameter a = O.OOSMeV^ 1 and values of the saturation energy per 
nucleon eo determined by a v and using the saturation density po = 0.1533/m -3 
for the standard nuclear matter. Rows (a) and (b) correspond to a v and a c values 
obtained from minimizing a 2 and x 2 /-/V respectively using both measured 2228+ex- 
trapolated 951 atomic mass excesses whereas rows (c) and (d) correspond to a v and 
a c values obtained from minimizing a 2 and "X 2 /N respectively using only the ex- 
perimentally measured 2228 atomic mass excesses. 





eo 
MeV 


a c 
MeV 


ro 
fm 


P 
fm 2 


C 


Kq 
MeV 


(a) 


15.260 


0.689 


1.254 


1.668 


2.07 


293.4 


(b) 


15.777 


0.710 


1.217 


1.676 


2.11 


301.1 


(c) 


15.409 


0.695 


1.244 


1.671 


2.08 


295.6 


(d) 


15.777 


0.710 


1.217 


1.676 


2.11 


301.1 



In table-3 the theoretical estimates of the pressure P and velocity of 
sound v s of standard nuclear matter have been listed as functions of nucle- 
onic density p and energy density e using the usual value of 0.005/MeV for 
the parameter a of energy dependence, given in eqn.(ll), of the zero range 
pseudo-potential. As for any other non-relativistic EOS, present EOS also 
suffers from superluminosity at very high densities. According to present 
calculations the velocity of sound becomes imaginary for p < 0.1/m, -3 and 
exceeds the velocity of light c at p > 5.3/Jo and the EOS obtained us- 
ing vn + TNI j!4j also resulted in sound velocity becoming imaginary at 
same nuclear density and superluminous at about the same nuclear density. 
But in contrast, the incompressibility Kq of infinite nuclear matter for the 
i>i4 + TNI was chosen to be 240 MeV while that obtained by the present 
theoretical estimate is 290-300 MeV which is in excellent agreement with the 
experimental value of Kq = 300 ±25 MeV obtained from the giant monopole 
resonance jll| and with the the recent experimental determination of Kq 
based upon the production of hard photons in heavy ion collisions which 
led to the experimental estimate of Kq = 290 ± 50 MeV • 

In Fig. [2J the plots of the velocity of sound v s in nuclear matter, the pres- 
sure P and the energy density e of nuclear matter as functions of nuleonic 
density p have been shown. The continuous line represents the velocity of 
sound in units of 10 _2 c for the standard nuclear matter, the dash-dotted 
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Fig. 1. The energy per nucleon e of nuclear matter as a function of p. The contin- 
uous line represents the curve for the present calculations using saturation energy 
per nucleon of -15.26 MeV and the dots represent the same using saturation energy 
per nucleon of -15.78 MeV whereas the dash-dotted line represents the same for the 
A18 model using variational chain summation (VCS) [13] for the spin and isospin 
symmetric infinite nuclear matter. 



line represents pressure in MeVfm~ 3 while the dotted line represents en- 
ergy density in MeV fm~ 3 for the standard nuclear matter. 

Keeping A constant while differentiating Eq.(4) and using Eq.(l) and 
setting the term dM nuc i eus (A, Z)/dZ \a equal to zero, one obtains 



Zstabl * = + WhereX = 2a >*» + K m »- m >W- (24) 

The second term in the numerator of the above equation is small compared 
to the atomic mass number A. The above relation which connects Z with A 
for the most stable nuclei can be written in a closed form as 



\A + 0.5aiA 2 / 3 } 

Z stable = : . 9/qi where a-y = a c /x, and a 2 = 4a sym /x. (25) 

[a 2 + aiA z l i \ 
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Table 3. Energy per nucleon e, pressure P, energy density e and velocity of sound 
v s as functions of nuclear density p using saturation energy per nucleon equal to 
— 15.260MeV and saturation density = 0.1533/m~ 3 for standard nuclear matter. 



p 


p/po 


e 


P 


e 








MeV 


MeVfm-' 3 


MeVfrn' 3 


in units of c 


.01 


.6523E-01 


-.7537E+00 


-.1677E-01 


.9382E+01 


.0000E+00 


.10 


.6523E+00 


-.1325E+02 


-.7633E+00 


.9257E+02 


.0000E+00 


.20 


.1305E+01 


-.1378E+02 


.2520E+01 


.1850E+03 


.2879E+00 


.30 


.1957E+01 


-.1138E+01 


.1689E+02 


.2813E+03 


.4700E+00 


.40 


.2609E+01 


.2341E+02 


.4829E+02 


.3849E+03 


.6207E+00 


.50 


.3262E+01 


.5896E+02 


.1020E+03 


.4989E+03 


.7442E+00 


.60 


.3914E+01 


.1048E+03 


.1830E+03 


.6262E+03 


.8443E+00 


.70 


.4566E+01 


.1605E+03 


.2958E+03 


.7696E+03 


.9248E+00 


.80 


.5219E+01 


.2254E+03 


.4447E+03 


.9315E+03 


.9895E+00 


.90 


.5871E+01 


.2993E+03 


.6340E+03 


.1114E+04 


.1042E+01 


1.00 


.6523E+01 


.3819E+03 


.8675E+03 


.1321E+04 


.1084E+01 



Table-4 provides values of ai, for different sets of co-efficients a c , a 
Table 4. The values of a\ and 0,2- 



a c 
MeV 


O'sym 

MeV 


a\ 


«2 


0.689 


22.209 


0.0153 


1.971 


0.710 


23.210 


0.0151 


1.973 


0.695 


22.435 


0.0153 


1.972 


0.710 


23.211 


0.0151 


1.973 



Since the values of a\ and 02 do not differ significantly for the different sets 
of co-efficient values a c and a sym , in Fig.|3]the plot of Z versus N=A-Z for the 
most stable nuclei is shown for the first set of values only. The continuous 
line represents the theoretical curve following the exact expression given by 
Eq.(25) while the dots represent the observed stable nuclei from the recent 
nuclide chart [T5| . 

The density dependence parameter (5 has the dimension of cross section. 
The density dependent term (1 — /3p 2 / 3 ) reduces the strength of the interac- 
tion which changes sign at high densities making it repulsive. It is a direct 
consequence of the Pauli blocking effect. Thus (1 — /3/3 2//3 ) can be interpreted 
as the probability of non-interaction arising due to the collision probability 
/3p 2 / 3 of a nucleon in nuclear medium of density p. The density dependence 
parameter (3 can be identified as the in medium effective nucleon-nucleon 
interaction cross-section Density dependence parameter (3 along with 
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0.0 0.2 0.4 0.6 0.8 1.0 

p (fnr 3 ) 

Fig. 2. The velocity of sound v s in nuclear matter, the pressure P and the energy 
density e of nuclear matter as a function of nuleonic desity p. The continuous 
line represents the velocity of sound in units of 10~ 2 c for the standard nuclear 
matter, the dash-dotted line represents pressure in MeVfm~ 3 while the dotted 
line represents energy density in MeVfm~ 3 for the standard nuclear matter. 



nucleonic density of infinite nuclear matter po can > therefore, provide the 
nuclear mean free path A = l/(po cr o)- Using value of the density dependence 
parameter (3 = 1.668/m 2 along with the nucleonic density of 0.1533/m -3 , 
the value obtained for the nuclear mean free path A is about 3.9/m. 



5. Summary and conclusion 

In summary, we conclude that five parameter least square fits to the 
Bethe-Weizsacker mass formula by minimizing a 2 and x 2 /N yield slightly 
different sets of energy coefficients a v , a s , a c , a sym and a p . Both a 2 and 
X 2 /N can not be minimized simultaneously by the identical sets of param- 
eters. The a 2 minimization is more appropriate since it reduces absolute 
errors involved in mass predictions. Identifying a v as the saturation energy 
per nucleon for the spin and the isospin symmetric cold infinite nuclear 
matter along with the value for saturation density equal to 0.1533/m~ 3 , 




50 100 150 200 250 

N 



Fig. 3. The plot of Z versus N for the most stable nuclei. The continuous line 
represents the theoretical curve while the dots represent observed stable nuclei. 

the nuclear incompressibility is found to be 290 — 300MeV which is in ex- 
cellent agreement with the experimental estimates from GMR as well 
as determination based upon the production of hard photons in heavy ion 
collisions |12| . The present theoretical estimate of nuclear incompressibility 
is in reasonably close agreement with other theoretical estimates obtained 
by INM ^Hj model, using the Seyler-Blanchard interaction or the rela- 
tivistic Brueckner-Hartree-Fock theory ^7]. The value of (3 provide nuclear 
mean free path A = 3.89 — 3.91/m which is in excellent agreement with the 
values derived by other methods [TS] , 
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